Ions at the air-water interface: An end to one hundred year old mystery? 
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Availability of highly reactive halogen ions at the surface of aerosols has tremendous implications 
for the atmospheric chemistry. Yet neither simulations, experiments, nor existing theories are able 
to provide a fully consistent description of the electrolyte-air interface. In this paper a new theory 
is proposed which allows us to explicitly calculate the ionic density profiles, the surface tension, and 
the electrostatic potential difference across the solution-air interface. Predictions of the theory are 
compared to experiments and are found to be in excellent agreement. The theory also sheds new 
light on one of the oldest puzzles of physical chemistry — the Hofmeister effect. 

PACS numbers: 61.20.Qg, 05.20.-y, 82.45.Gj 



Since van't Hoff's experimental measurements of os- 
motic pressure more than 120 years ago, electrolyte solu- 
tions have fascinated physicists, chemists, and biologists 
alike Q. The theory of Debye and Hiickel (DH) Q was 
able to address almost all of the properties of bulk elec- 
trolytes. On the other hand, electrolyte-air interface re- 
mains a puzzle up to now. The mystery appeared when 
Heydweiller Q measured the surface tension of various 
electrolyte solutions and observed that it was larger than 
the interfacial tension of pure water. While the depen- 
dence on the type of cation was weak, a strong variation 
of the excess surface tension was found with the type of 
anion. The sequence was reverse of the famous Hofmeis- 
ter series 0] , which was known to govern stability of pro- 
tein solutions against salting-out. An explanation for this 
behavior was advanced by Wagner [5| and Onsager and 
Samaras [|| (WOS), who argued that when ions approach 
the dielectric air- water interface, they see their image 
charge and are repelled from it. This produces a deple- 
tion zone which, with the help of thermodynamics, can be 
related to the excess surface tension. The theory and its 
future modifications Q, however, were unable to account 
for the Hofmeister series and showed strong deviations 
from the experimental measurements above lOOmM con- 
centrations. The fact that something was seriously wrong 
with the WOS approach was already clear in 1924, when 
Frumkin measured the potential difference across the air- 
water interface and found that for all halogen salts — ex- 
cept for fluoride — the electrostatic potential difference 
(air — water) was more negative for solution than for 
pure water [8(. This suggested that anions were able to 
approach the interface closer than the cations, or even be 
adsorbed to it! This contradicted the very foundation of 
the WOS theory. The confused state of affairs continued 
for the next 70 years, until the photoelectron emission ex- 
periments!^, [l(J, EH and the polarizable force fields sim- 
ulations [1 21 ] showed that Frumkin was right, and ions 
might be present at the interface. The situation, how- 



ever, remains far from resolved. Simulations predict so 
much adsorption that the excess surface tension of Nal 
solution becomes negative, contrary to experiments 
Furthermore, while the electron spectroscopy was finding 
the surface composition of solution to be enhanced in an- 
ions [icj |. vibrational sum-frequency spectroscopy (VSFS) 
indicated a significantly diminished anion population in 
the topmost layer [Tij ] . The two results appear to be con- 
tradictory. The questions, therefore, remain: Are there 
ions at the air-water interface? If so, why are they there 
and what are their concentrations? Besides its relation to 
the Hofmeister series, availability of highly reactive halo- 
gens at the surface of aerosol particles has a tremendous 
implication for the atmospheric chemistry and might help 
to explain the excessive rate of ozone depletion observed 
experimentally flB ]. In this Letter a theory will be pre- 
sented which allows all the pieces of this hundred year 
old puzzle to fit together. 

We begin by studying the excess surface tension 7 of 
an electrolyte solution. This can be calculated by inte- 
grating the Gibbs adsorption isotherm equation, d7 = 
— r+d/i+ — r_d/x_, where T± are the ion excess per unit 
area, and p± are the ionic chemical potentials. Suppose 
that the electrolyte is confined to a mesoscopic drop of 
water of radius R, corresponding to the position of the 
Gibbs dividing surface (GDS) [16(. Adsorptions are de- 
fined as 
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where p±(r) are the ionic density profiles and c& = 
p+(0) = P-(0) is the bulk concentration of electrolyte. 
If the system — water+vapor — contains N ion pairs, 
Eq. (HJ simplifies to T± = N/4irR 2 - c b R/3. 

Let us first consider the alkali metal cations such as 
Li + , Na + , and K + . These ions are small and strongly hy- 
drated. They can, therefore, be modeled as rigid spheres 
of hydrated radius ah and fixed charge q located at the 
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center. Water and air will be treated as uniform di- 
electrics of permittivities e w = 80 and e a = 1, respec- 
tively, with a discontinuity across the GDS. For a cation 
to move across the GDS, requires shedding its hydration 
sheath. For small, highly hydrated cations, this costs a 
lot of energy, resulting in a high potential energy barrier 
and a strong hardcore- like repulsion from the GDS. 

Suppose that a cation is located at position r p from 
the center of the drop. The electrostatic potential inside 
the electrolyte satisfies the usual DH equation, V 2 ip — 
K 2 tp = — Y^-S(r — r p ), where a — y / 8nq 2 Cb/e w kBT is the 
inverse Debye length. In the vapor phase the electrostatic 
potential satisfies the Laplace equation, V 2 tp = 0. For 
mesoscopic water drops of radius R ^> curvature can 
be neglected and the two partial differential equations 
can be solved using the Hankel transform [l7j]. Once 
the electrostatic potential is known, the work required to 
bring an ion from the bulk electrolyte to some distance 
z from the interface — z axis is oriented into the drop, 
with the GDS at z = — can be calculated using the 
Giintelberg charging process [lU . We find, 



W{z;a h ) = 
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2e w J s[s cosh(/ca/ l ) + k smh(fca/,,)] 



where s = \J k 1 + k 2 . Eq. ((2]) accounts for two fundamen- 
tal contributions: the interaction of an ion with its image 
across the interface, and for the loss of screening result- 
ing from breaking of spherical symmetry near the surface. 
In their theory of surface tension, Onsager and Samaras 
included ionic size by adopting the bulk spherically sym- 
metrical potential of Debye and Hiickel. By doing this, 
they have failed to account for the loss of screening near 
the interface, which leads to additional repulsion. This 
is one of the reasons why WOS theory significantly un- 
derestimates surface tensions of "hard" non-polarizable 
electrolytes such as NaF fl7j ]. 

While the alkali metal ions are strongly hydrated, the 
large halogen anions, such as iodine and bromide, have 
low electronic charge density and are weakly hydrated. 
To solvate an ion of radius ao , requires creation of a cav- 
ity and disturbance of the hydrogen bond network. For 
small cavities of radius ao < 4 A, the free energy cost 
scales with the volume of the void [l9(. If part of the 
ion leaves the aqueous environment, the cavitational en- 
ergy will decrease proportionately to the volume exposed. 
This results in a short range cavitation potential which 
forces ions to move across the air- water interface, 
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From bulk simulations |20[, we obtain v i=s 0.3fcsT/A 3 . 
For strongly hydrated alkali metal cations, the cavita- 
tional energy cost is too low to compensate the loss of 



hydration and the exposure of ionic charge to the low di- 
electric environment. For soft, unhydrated halogens, the 
situation is very different. As these ions move across the 
interface, they progressively shift their charge towards 
the part that remains hydrated, thus allowing them to 
reduce the cavitational energy at a small price in electro- 
static self energy 21 1 . 

To see how this works, consider a perfectly polariz- 
ablc ion modeled as a conducting spherical shell of ra- 
dius ao and charge g, free to distribute itself over its 
surface. The electrostatic self energy [2l| of such an ion 
when its center is located at distance z from the GDS 
is U s (z) — 2£ g - arccosu/io) — — ■ This expression is accu- 

rate for — ao < z < ao /4, and is exact for ions located pre- 
cisely at the GDS, z — 0. For such ions, U s (0) ~ q 2 /e w a Q , 
which is more than an order of magnitude lower than the 
electrostatic energy of a hard non-polarizable ion at the 
same position, ~ g 2 /4e a ao! 

Although fundamentally important at the interface, for 
distances z > ao, effects of ionic polarizability are neg- 
ligible. This can be verified by noting that for a hard 
ion located at z = ao, the electrostatic self energy is 
Uhard(ao) — 3q 2 /4e w ao. On the other hand, the self en- 
ergy of a perfectly polarizable ion at the same position is 
U so ft{a<o) = q 2 /2\n.{2)e w ao, which can be calculated ex- 
actly by resumming a series of images necessary to keep 
the ion at fixed potential. For ao ~ 2A, the difference be- 
tween Uhardio-o) and U so f t (ao) is about O.lfcsT. There- 
fore, for distances z > ag, the ionic polarizability can be 
neglected. The above calculation was performed in the 
infinite dilution limit. At finite ionic concentrations, po- 
larization effects will be even less significant, since all the 
induced interactions are doubly screened [22j . For z > ao 
the anion-interface electrostatic potential will, therefore, 
be well approximated by Eq. ([2]), with ah — * ao. The to- 
tal anion potential can then be obtained by interpolating 
between Eq. ([2]) and U s (z). We find 

W(z: a ) + val + for z > a 

Utot(z) = { w{a^ o)z + ^ + \ ay ( z) for o < z < a . (4) 
U s (z) + U cav (z) for - a < z < 

Since the electrostatic self energy is extremely large for 
z < — ao, no ion will be found at these distances. 

So far our discussion has been for perfectly polarizable 
ions. Real ions, however, have finite polarizability. The 
polarization potential for such ions has been derived in 
reference 2l|. For such ions, Eq. ^ should be modified 
by replacing the ideal potential U s (z), by the polarization 
potential U p (z) of reference 2l|. The potentials of all ions 
are plotted in Fig. [T] 

The ionic density profiles can now be calcu- 
lated numerically by solving the non-linear modi- 
fied Poisson-Boltzmann equation (mPB) V 2 0(r) = 



r[A] 



concentration [M] 



FIG. 1: The ion-surface interaction potentials at 1M salt con- 
centration. For hydrated ions there is a hardcore repulsion 
from the GDS located at r = 300A, and the cavitational en- 
ergy is omitted since it is constant. 



-^[p + {r)-p-(r% where 



p-(r) 
P+( r ) 



jR+a 47rr , 2 ^ el 3 q <t,(r)-(3U tot (r) 

N@(R - a h - r ) e -Pi<t>{r)-t3W{z;a h ) 



(5) 



r 

Jo 



R— at. 



47rr 2 dr e -Pi<f>( r )-f 3w ( z '' ah ) 



and O is the Heaviside step function. The excess surface 
tension of electrolyte solution can then be obtained by 
integrating the Gibbs adsorption isotherm equation (JTJ) 
with /3[i± = ln(c;,Aj_), where A± is the de Broglie thermal 
wavelength. 

We start with Nal. Since I - is large and soft, it should 
be unhydrated in the interfacial region. For its radius, 
we use the value calculated by Latimer, Pitzer, and Slan- 
sky [23| from fitting the experimentally measured free 
energy of hydration to the Born model. Latimer radii for 
halogens come out to be almost identical to the Pauling 
crystal radii, differing from them by only 0.1 A. Since 
in the bulk our theory reduces to the Born model, La- 
timer radii: aj — 2.26 A, as r — 2.05 A, aci = 1.91 
A, and ap — 1.46 A, are particularly appropriate. For 
ionic polarizabilities we will use the values from refer- 
ence Q: 7/ = 7.4 A 3 , ~/ Br = 5.07 A 3 , j C i = 3.77 A 3 , 
and jp = 1.31 A 3 . Our strategy will be to adjust the hy- 
drated radius of Na + to best fit the experimental surface 
tension for Nal. The same value of ah , will then be used 
to calculate the surface tension of other sodium salts and 
compare them to the experimental measurements [25j |. 
Fig. [2] shows the result of this procedure. We find that 
ah = 2.5 A gives an excellent fit to the experimental 
data for Nal. Since Br~ ion is also large and soft, we 
expect that it will also remain unhydrated in the interfa- 
cial region. This expectation is well justified, and a good 
agreement is obtained with the experimental data, Fig. 

The situation should be very different for F~, which is 
small, hard, and strongly hydrated. This means that just 



FIG. 2: Surface tension of NaF, NaCl, NaBr, and Nal. Na+ 
and Cl~ are partially hydrated with ah = 2.5 A and 2.0 A , 
respectively, F~ is fully hydrated with ah = 3.52 A, while the 
large halogens I - and Br - are unhydrated. Symbols are the 
experimental data from (2f|. 
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FIG. 3: Density profiles for NaF, NaBr, and Nal at 1M con- 
centration. The GDS is at r = 300A. 



like for cation, a hard core repulsion from the GDS must 
be explicitly included in the mPB equation. We then find 
an almost perfect agreement with the experimental data 
using the bulk hydrated radius of F~, ah — 3.52 A (2(|. 
The most difficult electrolyte to study theoretically is the 
usual table salt, NaCl. The size of chloride is sufficiently 
small that hydration must be taken into account. At the 
same time, it is also quite polarizable. We find that using 
Latimer size for Cl _ gives a very reasonable agreement 
with experiment, however, the agreement can be made 
perfect if we assume a small hydration, ah = 2.0 A. 

The density profiles for NaF, NaBr, and Nal are shown 
in Fig(3j In agreement with the polarizable force fields 
simulations, the density profiles for large halogens I - and 
Br~ differ significantly from the classical WOS picture. 
We find that there is a considerable concentration of an- 
ions at the GDS. However, unlike simulations [l3j, and in 
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agreement with the surface tension experiments, our ad- 
sorption always remains negative. The theory also agrees 
with the electron spectroscopy measurements showing 
that close to the GDS, there is a larger excess of anions 
over cations. However, just as was found using VSFS, the 
absolute concentration of anions at the surface is about 
half that of the bulk. We are, therefore, able to recon- 
cile the two sets of apparently contradictory experimen- 
tal results. Finally, we calculate the excess electrostatic 
potential difference across the interface, A%, for 1M so- 
lutions of NaF, NaCl, NaBr, and Nal . We obtain: +3.8, 
— 1.9, —9.1, and —14.0 mV, for the four salts respectively. 
These are quite close to the values originally measured 
by Frumkin [H, [2?|. In particular, one should note the 
change of sign of the electrostatic potential going from 
NaF to NaCl, first observed by Frumkin and confirmed 
by the present theory. 

We have seen that by adjusting only the hydrated ra- 
dius of sodium cation, we are able to account for the sur- 
face tensions of four different electrolyte solutions and for 
their values of A%. It should, therefore, be possible to 
predict the surface tensions and the Ax of other salts — 
as long as their anions are sufficiently large and weakly 
hydrated. This is the case, for example, for NaNC>3, 
NaIC>3, and NaClC>4. The only experimental data avail- 
able to us is for NaNC>3 [28[ which, once again, shows 
a good agreement between the theory and experiment. 
For NaC104, we find that at low concentrations the ex- 
cess surface tension is very small (slightly negative) but 
grows with increasing concentration of electrolyte. For 
1M NaC10 4 , the calculated value of A% is -42 mV, while 
the value originally measured by Frumkin was —48 mV 
27|. For 1M NaN0 3 and NaI0 3 , we obtain Ax = -8.2 
and —22 mV, respectively. 

We have presented a theory which allows us to cal- 
culate surface tensions and surface potentials for seven 
different electrolyte solution using only one adjustable 
parameter — the hydrated radius of sodium cation, ah = 
2.5 A. This value is very reasonable, lying between the 
Pauling crystal radius and the bulk hydration radius of 
Na + . In the case of five sodium salts for which there 
is experimental data available to us, the theory is found 
to be in very good agreement with the measured sur- 
face tensions. Using the same value of ah, we are also 
able to account for the experimentally measured electro- 
static potential differences across the solution-air inter- 
faces. The theory provides a very interesting picture of 
the interfacial region: alkali metal cations and fluoride 
anion are strongly hydrated and are repelled from the 
GDS. On the other hand, heavy halogens, Br - and I - , 
and the monovalent oxy-anions, NOg~, 10^ , and C107, 
are unhydrated, and as a result of their polarizability are 
significantly adsorbed to the surface. Nevertheless, their 
absolute concentration at the GDS remains below that 
of the bulk value. All these conclusions are in agreement 
with recent photoelectron spectroscopy and the VSFS 



measurements. In view of the success of the theory, it 
seems reasonable to hope that a fully quantitative under- 
standing of the Hofmeister effect might now be in sight. 

YL would like to thank C.J. Mundy for comments on 
the manuscript. This work is partially supported by 
CNPq, INCT-FCx, and by the US-AFOSR under the 
grant FA9550-06- 1-0345. 



[1] 
[2] 
[3] 



[5] 
[6] 

[7] 



[»] 
[9] 

[10] 

[11] 
[12] 



[13] 
[14] 
[15] 

[16] 

[17] 

[18] 
[19] 

[20] 

[21] 
[22] 

[23] 



Y. Levin, Rep. Prog. Phys. 65, 1577 (2002). 

P. W. Debye and E. Hiickel, Phys. Z. 24, 185 (1923). 

A. Heydweiller, Ann. Physik 33, 145 (1910). 

F. Hofmeister, Arch. Exp. Pathol. Pharmakol. 24, 247 
(1888); W. Kunz, P. Lo Nostra, and B. W. Ninham, Curr. 
Op. Colloid Interf. Sci. 9, 1 (2004). 

C. Wagner, Phys. Z. 25, 474 (1924). 

L. Onsager and N. N. T. Samaras, J. Chem. Phys. 2, 528 

(1934). 

A. L. Nichols III and L. R. Pratt, J. Chem. Phys. 80, 
6225 (1984). Y. Levin, J. Chem. Phys. 113, 9722 (2000); 
H. Ohshima and H. Matsubara, Colloid and Polymer Sci. 
282, 1044 (2004); D. S. Dean and R. R. Horgan, Phys. 
Rev. E 69, 061603 (2004); A. Onuki, J. Chem. Phys. 128, 
224704 (2008); M. Bier, J. Zwanikken, and R. van Roij, 
Phys. Rev. Lett. 101, 046104 (2008). 

A. Frumkin, Z. Physik. Chem. 109, 34 (1924). 

G. Markovich, S. Pollack, R. Giniger, and O. Chesh- 
novsky, J. Chem. Phys. 95, 9416 (1991). 

S. Ghosal and et al., Science 307, 563 (2005). 

B. Garrett, Science 303, 1146 (2004). 

L. Perera and M. L. Berkowitz, J. Chem. Phys. 95, 1954 
(1991); L. X. Dang and D. E. Smith, J. Chem. Phys. 99, 
6950 (1993); S. J. Stuart and B. J. Berne, J. Phys. Chem. 
A 103, 10300 (1999); P. Jungwirth and D. J. Tobias, J. 
Phys. Chem. B 106, 6361 (2002); P. Jungwirth and D. J. 
Tobias, Chem. Rev. 106, 1259 (2006); D. Horinek and 
R. R. Netz, Phys. Rev. Lett. 99, 226104 (2007); M. A. 
Brown and et al., Phys. Chem. Chem. Phys. 10, 4778 
(2008). 

T. Ishiyama and A. Morita, J. Phys. Chem. C 111, 721 
(2007). 

E. A. Raymond and G. L. Richmond, J. Phys. Chem. B 
108, 5051 (2004). 

E. M. Knipping, M. J. Lakin, K. L. Foster, P. Jung- 
wirth, D. J. Tobias, R. B. Gerber, D. Dabdub, and B. J. 
Finlayson-Pitts, Science 288, 301 (2000). 

C. -H. Ho, H.-K. Tsao, and Y.-J. Sheng, J. Chem. Phys. 
119, 2369 (2003). 

Y. Levin and J. E. Flores-Mena, Europhys. Lett. 56, 187 
(2001). 

E. Z. Guntelberg, Z. Phys. Chem. 123, 199 (1926). 

K. Lum, D. Chandler, and J. D. Weeks, J. Phys. Chem. 

B 103, 4570 (1999). 

S. Rajamani, T. M. Truskett, and S. Garde, PNAS 102, 
9475 (2005). 

Y. Levin, Phys. Rev. Lett 102, 147803 (2009). 

M. E. Fisher, Y. Levin, and X.-J. Li, J. Chem. Phys. 101, 

2273 (1994). 

W. M. Latimer, K. S. Pitzer, and C. M. Slansky, J. Chem. 
Phys. 7, 108 (1938). 



■5 



[24] N. C. Pyper, C. G. Pike, and P. P. Edwards, Molecular 
Phys. 76, 353 (1992). 

[25] N. Matubayasi, K. Tsunemoto, I. Sato, R. Akizuki, 
T. Morishita, A. Matuzawa, and Y. Natsukari, J. Col- 
loid and Interface Sci. 243, 444 (2001). 

[26] E. R. Nightingale Jr., J. Phys. Chem. 63, 1381 (1959). 



[27] J. E. B. Randies, Advances in Electrochemistry and Elec- 
trochemical Engineering (Interscience, New York, 1963), 
vol. 3, p. 1. 

[28] N. Matubayasi and R. Yoshikawa, J. Colloid and Inter. 
Sci. 315, 597 (2007). 



